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Abstract 

The multiplicity distribution of particles in relativistic gases is studied 
in terms of Tsallis' nonextensive statistics. For an entropic index q > 1 
the multiplicity distribution is wider than the Poisson distribution with 
the same average number of particles, being similar to the negative bino- 
mial distribution commonly used in phenomenological analysis of hadron 
production in high-energy collisions. 

1 Introduction 

It is always amusing to recognize that the concepts of equilibrium statistical 
physics, and hence of thermodynamics, can be applied in (at first sight) unex- 
pected branches of physics. For example, it has long been known that several 
features of multiparticle production in high-energy hadronic collisions are well 
described by thermostatistical models JjJ ^, [jj . More recently, it has been shown 
that a simple thermal description of particle abundances works quite well for 
the hadrons produced in elementary electron-positron collisions (J, ||. This 
indicates that the hadronic system reaches statistical equilibrium, at least for 
some observables. One might try to justify this in terms of the large number of 
degrees of freedom in the final state. 

Thermostatistics is also an important issue in relativistic heavy ion physics. 
The main focus of recent investigations in this field is the production and iden- 
tification of the quark-gluon plasma, a phase of nuclear matter consisting of 
unconfined quarks and gluons. Thermostatistical aspects play a central role in 
these investigations, as it is expect that some form of statistical equilibrium is 
attained in high-energy collisions of two heavy nuclei, at least in the central 
region §,0. 

One of the signatures of "thermal" multiparticle production is the exponen- 
tial form of the transverse energy distribution of the produced hadrons. We are 
tempted to interpret the slope parameter of this distribution in terms of a tem- 
perature of the final state. In fact, up to the ISR energy region, the transverse 
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energy distribution is consistent with such an interpretation. The same happens 
to the observed hadron abundances. 

However, at higher energies, the thermal interpretation of the transverse 
energy spectra should be modified. The SPS experiments revealed a significant 
deviation from the exponential transverse energy distribution, together with 
the violation of several scaling laws. Usually, these aspects are understood 
in terms of the onset of hard QCD processes such as the formation of mini- 
jets. Furthermore, even if local thermal equilibrium is attained, dynamical 
effects such as collective flow may entangle with the thermal transverse spectra. 
Consequently, for these energies, a pure statistical description to the overall 
system cannot be applied. 

Recently, it was pointed out H, ^|, [l0| that the nonexponential behav- 
ior of the transverse energy distribution in high-energy multiparticle production 
processes (including electron-positron annihilation) can be described rather well 
by the nonextensive thermostatistics of Tsallis [fijl. In addition to the excel- 
lent quality of the fit to transverse energy spectra, an interesting point of the 
extended statistical approach to multiparticle production is that the "temper- 
ature" stays almost constant, independent of the incident energy, conveniently 
recovering the classical Hagedorn scenario [0, 0, [§]. This opens the possibil- 
ity of reviving the statistical description of the complex dynamics of hadronic 
collisions in a unified way, including the onset of semi-hard QCD processes, pro- 
vided that we generalize the concept of statistical equilibrium. Differently from 
the standard statistical mechanics of Boltzmann and Gibbs, the Tsallis nonex- 
tensive thermostatistics contains one extra parameter, the "entropic index" q, 
originated in the definition of the entropy. The physical origin of this parameter 
is not yet completely understood and a lot of discussions are in course. It can be 
attributed to various factors, like the dynamical suppression of certain domain 
of the phase space (multifractal structure near the critical point), the presence 
of long range forces, or fluctuations in a small system |l3|, Oj. For an entropic 
index q = 1 Tsallis statistics reduces to the Boltzmann-Gibbs theory If q > 1, 
rare events are enhanced relative to the Boltzmann-Gibbs case. If g < 1, fre- 
quent events are privileged. So far, all phenomenological applications of Tsallis 
statistics to high-energy collisions have found q > 1 in these processes. Some 
discussion on the possible origin of nonextensive behavior in hadronic systems 
can be found in Refs. [jl5[ ^6| . 

The transverse energy distribution is not the only probe of the formation 
of a locally thermalized source of particles in high-energy collisions. Another 
important observable is the multiplicity distribution of the produced hadrons. 
The negative binomial distribution Q 

T(N + k) (N/k) N 

N r(jv + i)r(fc) {i + N/k) N + k ' 1 ; 

is often used to express hadron abundances in high-energy processes, where N 
is the number of particles and N the average multiplicity. The k parameter is 
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related to the variance D 2 = iV 2 — N of the distribution by 

n — N 2 

D 2 = N+—. (2) 

Although the negative binomial distribution is known to arise from some specific 
processes, such as Bose- Einstein particles with different sources |l|], its use in 
high-energy multiparticle production is rather phenomenological. 

In this paper, we discuss how the generalized statistical mechanics of Tsallis 
(with q > 1) affects the multiplicity distribution in hadronic collisions. We show 
that for q close to unity the distribution of particle numbers in a high-energy gas 
is approximately negative binomial. The Boltzmann-Gibbs and Tsallis statistics 
are applied to a simple meson gas model, and the results are compared to 
experimental data. 



2 Nonextensive q statistics 

First, let us review briefly Tsallis' statistics and derive some formulas necessary 
for our discussion. Tsallis generalized the usual Boltzmann-Gibbs thermostatis- 
tics by introducing the "g-entropy" 

l-Y v q 

S = _ i , (3) 
where p a is the probability associated with microstate a, normalized as 

X> = 1 - (4) 

a 

Tsallis also introduced "g-biased" averages of observables, 

(0} = i^O oP «, (5) 

y a 

the normalization factor being 

C q = J2Pa- ( 6 ) 

a 

The g-biased microstate probability 

Pa = £ (7) 
Q 

is the probability to be used in the calculation of physical quantities. The 
entropic index q is a real number, and in the limit q — * 1 the Boltzmann-Gibbs- 
Shannon entropy is recovered. In this limit the g-average also reduces to the 
usual one. The equilibrium probabilities p a are determined by maximizing the 
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entropy under appropriate constraints. These are the normalization condition 
([|), and the fixed (average) value of the energy E and of any other relevant con- 
served quantity. We assume for simplicity that there is only one such quantity, 
a conserved charge Q. The variational principle is then 

SS + a ^2 s Pa ~ Pt SE + 7 6Q = , (8) 

a 

where 

E = J2EaPa, (9) 

a 

Q = Y,Q<3«> ( 10 ) 

a 

and E a and Q a are the energy and charge of the microstate a. The constants 
a, [3t and 7 are Lagrange multipliers. The later two are associated with the 
temperature T and chemical potential fi (in the sense of the second law of 
thermodynamics) as 

where V is the volume occupied by the system. Solving the variational equation 
(||) we obtain the Tsallis distribution, 

p a = ±-{exp 9 [-0(E a -nQ a )]} q , (13) 

where we have defined the "g-exponential function", exp 9 , by 

exp^^l-te-lM]- 1 /^ 1 ). (14) 
We have also introduced the generalized partition function Z q as 

Z g ((3,(x,V) = ^2{vq> q [-0(E a -nQ a ))} q . (15) 

a 

Note that in the limit of q — > 1 we have 

lim exp (A) = e A (16) 

so that Eqs. jl^ ) and ( p!5| ) reduce to the corresponding quantities of the usual 
Boltzmann-Gibbs statistical mechanics. The parameter [3 appearing in these 
equations is not the inverse temperature Lagrange multiplier f3r — l/T; (3 is 
related to temperature T by 

T= /3- 1 + (q-l)(E-fiQ) 
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The Tsallis non-intensive temperature T meets some difficulties when confronted 
with the zeroth law of thermodynamics. This point has been investigated re- 
cently Jl9| , and it was shown that the quantity 

T = [3- 1 + (q-l)(E-[iQ) , (18) 

sometimes called the "physical" temperature, can provide a better characteri- 
zation of thermal equilibrium. 

Instead of using the chemical potential fx to control the average charge Q, as 
in the grand canonical approach outlined above, it is sometimes necessary to im- 
pose charge conservation exactly. This is particularly important for small values 
of Q, when fluctuations around the average become significant. Then a canoni- 
cal treatment is preferable, in which case the generalized partition function for 
a fixed charge Q is given by 

Z g {p, Q, V) = $ (Q - Qa) [ex P? (~(3E a )] q (19) 

a 

where S(Q — Q a ) is a Kronecker delta. The generalized canonical probability is 

[•*.<-«■>]•■ <*» 

3 Integral representation of the Tsallis distribu- 
tion 

For further development of the theory, it is useful to introduce the following 
integral representation, valid for q > 1, 

/>oo 

[exp ? (A)] 9 =/ dxG(x)e xA , (21) 
Jo 

where 

G[x) = T§) e ' vx (22) 

with v = l/{q-l). Since G (x) > and 

G(x)dx = l, (23) 



the function G (x) can be considered as the probability distribution of the vari- 
able x £ [0, oo). The maximum of G(x) is at x = 1, and the moments of x 
are 

" x-G (x) dx = (" + ")fr + "-l)-fr + D . (24) 

v n 

In particular, we have 

x=^±± = q, (25) 
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x 1 = — 2 = q(2q-l), (26) 

so that the width of G{x) is 

a x = \j^-x 2 = Vg(g-i) • (27) 

For q — > 1, G (x) tends to the Dirac delta function 5 (x — 1). 

Using the integral representation (|2l]) , we can express the generalized parti- 
tion function (pUf ) as 



Z q {fl,n,V)= / dxG(x)Z(x0,n,V), (28) 
Jo 

where 

Z(/5,M ) y)=^exp[-/5(£ (l -A)] (29) 

a 

is the Boltzmann-Gibbs partition function. 

The integral representation of the q-exponential function is useful for devel- 
oping approximations when q is close to unity. As we have seen, in this case the 
function G (x) is sharply peaked at x = 1. Thus, in any expression of the form 

dxG (x) e^ x) , (30) 

we may expand the exponent / (x) around x = 1 as 

/ (a) ~/(l) + /'(l) Or - 1) + ••• , (31) 

and obtain 

j „ e /(D-/'(D / dxG (x)e xf 'W 

= e /w-m { ex Pg [/'(i)]} 9 

~ exp{/(l) + (g-l)/'(l)[l + /'(l)/2]} • (32) 

In the last step we have used 

[ex Pq (A)] q ^ex V [A+(q-l)(A + A 2 /2)} , (33) 

valid for q ~ 1. Applying these approximations to Eq. (p8|), the generalized 
partition function Z q {(3, [a, V) can be written in terms of the Boltzmann-Gibbs 
function Z{(3, fi, V) as 

Z q (p, fi, V) ~ Z(/3, (j,, V) exp {( ? - 1) [-P(E - fiQ) + [3 2 {E - [iQ) 2 /2] } , 

(34) 

where 

<91nZ 

E-(iQ = ——. (35) 

Integral representations of exp 9 also exist for q < 1 (2fj, |2l| . We will not 
explore these here because, as mentioned in the Introduction, q > 1 seems to 
be the case of interest in high-energy collisions. 
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4 Particle multiplicity distribution 



The probability that the system has exactly N particles is given, in the Tsallis 
statistics, by 

P N = 5{N - N a ) p a (36) 

a 

where N a is the particle number in state a. Defining the partition function 
restricted to N particles, 

I 9 



Z\ N \(3, n,V)=Y, S(N - N a ) {exp 9 [-0 (E a - nQa)]}' 



(37) 



we have 



Pn = 



z„ 



(38) 



Using the integral representation of the ^-exponential function, we can write 
zW(fl,iM,V)= dxG{x)ZW (x/3,n,V) 



(39) 
Jo 

where Z^ N ^ is the Boltzmann-Gibbs restricted partition function 

03, i x,V)=Y J S(N ~ N a ) exp [-0 (E a ~ f,Q a )] . (40) 

a 

It is useful to introduce the generating function for the multiplicity distri- 
bution, defined by 



F(t) = £;t"p w 



AT=0 

For example, if Pn is a Poisson distribution, 

then its generating function is an exponential, 

F P (t) =exppV(i- 1)]. 
For the negative binomial distribution, Eq. (pi), we get 



F NB (t) = 



= ex Pg [N (t - 1)] 



(41) 

(42) 
(43) 

(44) 



where q — 1 + 1/k. It is suggestive to note that the negative binomial generating 
function is obtained substituting the exponential in the Poisson generating func- 
tion by the g-exponential. For large values of k the negative binomial generating 
function has the asymptotic behavior 

-2 



F NB (t) ~ exp 



(45) 



a result that will be useful later on. We see that in the limit k — ► oo the negative 
binomial reduces to the Poisson distribution. 
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5 Relativist ic ideal gas 



In order to study the Tsallis statistics of a hadronic gas, it is tempting to apply 
the results of the previous sections to an ideal relativistic Boltzmann gas. Let 
us consider h different particle species, with masses m, and charges qi, i = l...h. 
In this case, the Boltzmann-Gibbs partition function is 

Z M ,y)=exp^^ exp^iiqA , (46) 

where 

Above, gi is the statistical factor of particle i and K2 (z) is a modified Bessel 
function. 

The iV-particle partition function is 

/ h \ N 

ZW{0,IM,V) = ^ (vJ2HP) exp(/?M%) J (48) 



i=i 



and, from Eq. (|38|), the multiplicity distribution of the ideal Boltzmann gas is 
seen to be Poissonian. The average number of particles is 

N = n(p,n)V, (49) 

where 

h 

n((3, fi)=Y / MP) exp(Au%) (50) 

i=l 

is the particle density. 

We meet a problem when we try to apply Tsallis' statistics to the ideal 
gas: for q > 1 the generalized partition function Z q is infinite. The integral 
representation of Z q , Eq. (p8|), diverges if the ideal gas partition function ( ff6| ) 
is substituted in the integrand. This happens because for (3 — > 0, 

**C0 - ^3 . (51) 

so that Z(x(3, /j,, V) has an essential singularity at x = which cannot be re- 
moved by any power of x in G (x). Therefore, the g-statistics of an ideal rela- 
tivistic gas cannot be defined — there is no ideal Tsallis gas with q > 1. 



6 Relativistic van der Waals gas 

When restricted to N particles, the generalized ideal gas partition function 
Zq N ^ has a well defined integral representation in terms of the corresponding 



8 



Boltzmann-Gibbs function (Eq. fl4g|)), provided 

3(g — 1) 

It is the production of particles above this limit that causes the divergence in the 
partition function of the Tsallis ideal gas. High multiplicities can be suppressed 
if a repulsive interaction exists among the produced particles. A simple way 
to include such interactions is to introduce a "van der Waals" excluded volume 
simulating the effect of hard-core potentials among particles. The van der Waals 
gas model has been frequently used to analyze particle abundances in high- 
energy heavy ion collisions p2l 

Let vq be the excluded volume associated to a particle. The corresponding 
partition function for N particles is obtained by replacing the volume of the 
system V by V — Nv , 

Z<"> (/?, n, V) -> (f3, M , V - Nv ) 8 (V - Nv ) . (53) 

The Heaviside 0-function limits the number of particles inside the volume V to 
N < V/vq. The partition function for the relativistic van der Waals gas is then 

Z (/?, H,V) = J2 ») N (V ~ Nvof 8 (V— Nvo) (54) 

N 

where n(/3, fx) is the ideal gas density given in Eq. (|50"|). In the large V asymptotic 
limit, this can be written as 

Z ((3, n, V) = exp I ^-W [v n((3, ft)] J , (55) 

where W(x) is the Lambert function [^3[ defined by the equation 

W{x)e w{x) = x . (56) 
Similarly, the generating function of the excluded volume gas is found to be 

Z{p,fi,V) [v J 
The Lambert function W (z) has the asymptotic behavior, in the limit z — » 

oo, 

W (z) ~ Inz + lnlnzH . (58) 

The essential singularity at (3 — ► of the ideal Boltzmann gas partition function 
is thus reduced to a power one, 

Z{p^,V)^r Who ■ (59) 

Therefore, the integral in 

Z q (i3,Li,V) = J dxG(x)expl^-W[v n(xf3,Li)]\ , (60) 
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converges as long as 

*<1 + ^, (61) 

as expected from Eq. (]52|), since the hard core suppresses the production of 
particles above the maximum number V/vq. In this range of q values, the 
generating function of the multiplicity distribution in the Tsallis - van der Waals 
gas is given by 

F(t) = - 1 / dxG(x)exp\-W[tv n(x(3, f i)}\ . (62) 

Zq {P, H, V) J [V J 

7 Tsallis and van der Waals corrections to the 
ideal gas 

In order to examine in more detail the effect of Tsallis statistics on the multi- 
plicity distribution, let us consider the case in which q — 1 and vq are both small 
(respecting the limit of Eq. (|6l|)). The Lambert function W(z) has the series 
expansion p3[ 

00 ( \n—l 

^) = E iZ ^T— ( 63 ) 

n=l 

so that for — > we have 

Vy(two^) ^ tv n - (to n) 2 H (64) 

and the generating function for the van der Waals-Tsallis relativistic gas can be 
written as 

1 r°° 

F{t)~— I dxG (x) exp {tVn(x/3, fj) [1 - tv n{x(3 7 fi)] } . (65) 

Zq Jo 

For q — ► 1, this integral can be calculated with the help of Eq. j32|) . To first 
order in q — 1 and i>o we obtain 

F{t) ~ exp{(* - l)Vn[l + (g - l)A(7nA - 1) - 2v n] 

+ (t - l) 2 (Vn) 2 [(g - l)A 2 /2 - vo/V}} (66) 

where 

A(M = -^- (67) 

Comparison to Eq. (^5|) shows that this is the generating function of a negative 
binomial distribution, with average number of particles 

77 = Vn [l + (q- l)\(Vn\ - 1) - 2v n] (68) 

and a (large) fc-parameter given by 

\ = {q-l)\ 2 -2 V f. (69) 
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Strictly speaking, a negative binomial distribution is obtained only if k > 0, or 



•> 1 + w- <*» 

In most cases of interest for hadron production the value of A(/3,/i) is large, so 
that Eqs. (]70|) and @ can both hold. If q is bellow the limit of Eq. @, fc is 



negative, and the multiplicity distribution will have a binomial-like form. From 
Eq. (|6^) we see that the effect of Tsallis statistics is to make the multiplicity 
distribution wider than the Poisson distribution with the same average number 
of particles. The effect of the van der Waals excluded volume is to reduce the 
width relative to Poisson. 



8 Charged particle multiplicity in hadronic col- 
lisions 

Pions are the most copiously produced particles in high-energy hadronic colli- 
sions. They come not only from the thermally equilibrated gas, but also from 
the decay of heavier hadrons formed in the gas. A precise description of par- 
ticle multiplicities has to take into account how much these unstable hadrons, 
called resonances, participate in the production of the observed particles. In the 
present analysis, we intend to understand the basic influence of the g-statistics 
on the particle multiplicity distribution. Therefore, to simplify the picture, we 
will study the behavior of the lowest-lying non-strange mesons, it, rj, p and lu in 
thermal equilibrium. The 77, p and lo mesons are resonances which decay into 
pions. The essential features of the baryon-free hadronic gas can be described 
by the mixture of these mesons, and the major conclusions of the present work 
will not be changed by this simplification. 

In Fig. [l], we show the multiplicity distribution of charged particles pro- 
duced in pp collisions at momentum pi a b — 250 GeV/c. Triangles are data 
from the NA22 experiment psj] . The average number of charged particles is 
N = 7.88 ± 0.09 and the width of the distribution is D = 4.10 ± 0.05. The lines 
are multiplicity distributions calculated for the meson system described above, 
treated as an ideal Boltzmann gas. Two different temperatures were considered 

T = 190 MeV (solid line), and T = 160 MeV (dashed line) — correspond- 
ing to the range of values found in analyses of pp reactions [^6| . In both cases 
the volume of the gas was chosen such as to reproduce the measured average 
multiplicity. For the higher temperature this gives V = 18.8 fm 3 , and for the 
lower, V = 43.0 fm 3 . We assumed that the meson gas has a total electric charge 
Q = 1, which is the typical value in pp collisions. Changing this charge to or 
2 has only a small effect on the multiplicity distribution. All calculations were 
performed in the canonical framework. 

Due to resonance decay, and exact charge conservation in the canonical en- 
semble, the calculated multiplicity does not follow exactly the Poisson distribu- 
tion obtained in the grand canonical approach of Sec. o. The deviation from 
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Poisson is conveniently measured by the quantity 



- = — - (71) 

a definition motivated by the /c-parameter of the negative binomial distribution 
(compare with Eq. (||)). The ideal gas distributions shown in Fig. [t]have k ~ 40 
and are, therefore, wider than Poisson. They are well approximated by a neg- 
ative binomial distribution with the same k parameter. The large value of k 
shows that the deviation from Poisson caused by resonance decay and charge 
conservation is relatively small. This is reflected in Fig. |l|, where multiplicity 
distributions constrained to have the same average present similar widths, de- 
spite the very different gas temperatures and volumes. It also indicates that 
the ideal gas model cannot explain the multiplicity data, in particular its large 
width. For all reasonable values of T and V consistent with the experimental 
average multiplicity, the ideal gas shows approximately the same narrow distri- 
bution. One may think that including more hadronic resonances in the model 
improves the situation. However, Becattini et al. || find k > 10 in their analy- 
sis of e + e~ data, which includes many resonances. Hadronic collisions produce 
somewhat wider multiplicity distributions; the NA22 data shown in Fig. [TJ have 
k = 7, and as the energy increases k becomes smaller p7| . We thus conclude 
that the multiplicity distribution of the ideal relativistic gas is not consistent 
with the hadron-hadron experimental data. 

The situation doesn't get better if the meson system is treated as a rel- 
ativistic van der Waals gas. This is seen in Fig. |^, where the NA22 data is 
compared to the canonical multiplicity distribution of the meson gas with ex- 
cluded volume v = 0.368 fm 3 . Again two temperatures, T = 190 MeV (solid 
line) and 160 MeV (dashed line), were considered, with volumes V = 24.5 fm 3 
and 49.2 fm 3 , respectively, determined by the experimental average multiplicity. 
The excluded volume reduces the width of the multiplicity distribution, relative 
to the ideal gas with the same average multiplicity; we have k ~ 100 for the 
van der Waals gas distributions shown in Fig. |^. The width reduction occurs 
for any value of vq. Therefore, the van der Waals gas model cannot explain the 
experimental multiplicity distribution; it is even less satisfactory than the ideal 
gas model. 

Tsallis statistics for q > 1 produces multiplicity distributions that are wider 
than the Boltzmann-Gibbs ones (see Sec. 0). In Fig. [|, we show some examples 
of the canonical multiplicity distribution obtained in g-statistics. The solid 
line corresponds to a gas with = 190 MeV, V = 9.13 fm 3 , q = 1.0097, 
and the dashed line to /T 1 = 160 MeV, V = 31.0 fm 3 , q = 1.0030. In both 
calculations the excluded volume is vq = 0.368 fm 3 . We note that a even a small 
deviation of q from unity causes a substantial broadening of the distribution, 
taking it much closer to the experimental data. Thus, it seems possible to give 
a thermostatistical description of high-energy hadronic multiplicities, provided 
nonextensivity effects are taken into account. 



12 



9 Concluding remarks 



In the present work we have studied the effect of the Tsallis nonextensive statis- 
tics (with q > 1) on the multiplicity distribution of particles in a high-energy 
gas. In order to avoid divergences in the generalized partition function of the 
ideal gas, a van der Waals hard-core interaction was introduced as a physical 
regularization procedure. Nonextensive statistics gives rise to multiplicity dis- 
tributions that are broader than the Boltzmann-Gibbs ones. In the small q — 1 
limit, the Tsallis multiplicity distribution can be approximated by a negative 
binomial distribution, with fc-parameter proportional to (q — 1) . We devel- 
oped a simple meson gas model in order to compare the different statistics to 
experimental data from hadronic collisions. The multiplicity distribution ob- 
tained with the Boltzmann-Gibbs statistics has a too narrow width compared 
to the experimental results. A small increase of q from unity makes the distribu- 
tion considerably wider, and a good agreement with the data can be obtained 
within Tsallis statistics. These results provide an interesting perspective on 
why hadronic multiplicities have a negative binomial distribution, and suggest 
that nonextensive thermal phenomena play an important role in high-energy 
collisions. 

In order to get a more complete picture of the nonextensive thermostatistical 
effects in hadron production, we need to extend our analysis using a larger 
resonance gas, and including other observables such as particle ratios, transverse 
momentum spectra, and correlations. Work along these lines is in progress. 

We wish to thank C. Tsallis and G. Wilk for fruitful discussions. This work 
was partially supported by FAPERJ and CNPq. 
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Figure 1: Multiplicity distribution of charged particles produced in pp collisions 
at 250 GeV/c. The lines are calculations based on the ideal Boltzmann gas 
model. The solid line corresponds to T = 190 MeV and V = 18.84 fm 3 , the 
dashed line to T = 160 MeV and V = 42.96 fm 3 
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Figure 2: Charged particle multiplicity distribution in pp collisions at 
250 GeV/c. The lines are multiplicity distribution of van der Waals gases with 
T = 190 MeV, V = 24.46 fm 3 (solid), and T = 160 MeV, V = 49.21 fm 3 
(dashed). In both calculations the excluded volume is vq = 0.368 fm 3 . 
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Figure 3: Charged particle multiplicity distribution in pp collisions at 
250 GeV/c. The lines are multiplicity distribution of Tsallis gases with 
/T 1 = 190 MeV, V = 9.13 fm 3 , q = 1.0097 (solid), and /3' 1 = 160 MeV, 
V = 31.0 fm 3 , q = 1.0030 (dashed). In both cases v = 0.368 fm 3 . 
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